Double-well magnetic trap for Bose-Einstein condensates 
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We present a magnetic trapping scheme for neutral atoms based on a hybrid of loffe-Pritchard and 
Time-averaged Orbiting Potential traps. The resulting double-well magnetic potential has readily 
controllable barrier height and well separation. This offers a new tool for studying the behavior of 
Bose condensates in double-well potentials, including atom interferometry and Josephson tunneling. 
We formulate a description for the potential of this magnetic trap and discuss practical issues such 
as loading with atoms, evaporative cooling and manipulating the potential. 
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Bose-Einstein condensates (BECs) of dilute alkali gases 
have been the subject of a great deal of attention since 
they were first realized ||, ||, ^ . One area in which there 
has been considerable theoretical interest is the behavior 
of Bose condensates confined in double-well potentials. 
The process of splitting a condensate in two has been in- 
vestigated by raising a potential barrier in the center of a 
harmonic trap |Q |5[ |^, |^, ^ . This work has given insight 
into the dynamic behavior of condensates and modifica- 
tion of collective excitations. Also extensively studied is 
the area of Josephson junctions |lO[ |l^, |l^, p^ , 
where the condensate mean field affects the tunneling 
between wells and macroscopic quantum self-trapping is 
possible. Other researchers have proposed experiments 
involving the relative phase of the two condensates. For 
example, uniting the relative phase of two initially inde- 
pendent Bose condensates [|l^, investigating quan- 
tum and thermal fluctuations of the phase and de- 
tection of weak forces Q . 

To date the only experimental demonstration of a 
double-well potential for Bose condensates has been a 
hybrid configuration involving a harmonic magnetic po- 
tential divided into two regions by a repulsive optical 
field [|l9j. This scheme was used to observe interference 
fringes of two overlapping Bose condensates, in the first 
demonstration of long range coherence. More recently, 
the same configuration was used to measure the presence 
of vortices induced in one of the wells . One feature 
of the apparatus is that there is high degree of stability 
between the magnetic and optical fields. 

In this work we address the problem of producing a 
purely magnetic double-well potential suitable for study- 
ing BECs. Our scheme is an extension of the existing 
loffe-Pritchard (IP) and Time-averaged Orbiting Poten- 
tial (TOP) traps used in many BEG experiments. We 
begin by briefly summarizing the IP trap and review how 
a double quadrupole-like potential can be formed. Bose 
condensates cannot be confined in this potential, because 
of Majorana spin-flip loss, but we show how the addition 
of a rotating bias field can resolve this problem. The 
focus of this work is the development of a theoretical de- 
scription of the resulting double- well magnetic potential. 



Practical issues such as loading, evaporative cooling and 
controlling the shape of the potential are considered in 
detail. 



I. lOFFE-PRITCHARD TRAPS 

loffe-Pritchard (IP) traps |2l|, ^ have been used ex- 
tensively in the realization of BEG in alkali gases. The 
basic configuration for an IP trap is illustrated in Fig. [|. 
Four long "bars" with currents in alternate directions 
run parallel to the z axis. These generate a quadrupole 
field in the x-y plane with gradient B' , which leads to ra- 
dial confinement. Two "pinch" coils have currents in the 
same direction and are spaced to give a harmonic local 
minimum of the axial magnetic field with curvature B" . 
These coils provide axial confinement but also a large 
bias field Bp. Finally there are the two "nulling" coils 
(usually in the Helmholtz configuration) which provide a 
uniform axial field i?„. These allow the total bias field 




FIG. 1: The layout of coils and bars for a loffe-Pritchard mag- 
netic trap with current direction indicated by arrows. The 
four bars lie on corners of a square. Current in each bar flows 
in the direction opposite to that of its closest neighbors, and 
the magnitude of the current is the same for each bar. The 
two pairs of coils with large and small radii are the nulling 
and pinch coils respectively. The currents flow in the opposite 
sense so that they produce opposing magnetic fields. 
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TABLE I: A summary of reported IP trap parameters. In 
many cases the traditional geometry of Fig. m is replaced by 
a variant with some specific advantages. References denoted 
* have field quantities calculated from the trap oscillation 
frequencies, and f indicates an equivalent radial gradient for 
traps that are harmonic. 



at the trap center, 



Bn — Bn 



Br, 



(1) 



to be reduced. For BEC experiments this is important 
because it leads to ti ght er radial confinement and simph- 
fies RF evaporation |23| . 

The magnetic field components for this geometry, to 
second order about the trap center, are given by the equa- 
tions 



Bx — 

By = 

B. = 



+B'x -^B"xz, 
-B'y -W'yz, 
Bo +\b"{z^ 



1 ' n 



(2) 



where r is the radial coordinate (r^ = x^+y^). The terms 
with coefficient B' are generated by the four bars. Atoms 
in a weak-field seeking state are trapped at the minimum 
of the field magnitude. For small thermal atomic clouds 
and Bose-Einstein condensates (fc^T < iiBo), the field 
magnitude is calculated by a binomial expansion to sec- 
ond order, to give 



1 B" 



-B"z\ 



(3) 



where Bo > 0. The radial curvature is large when 
B'^ /Bo ^ B" /2, which is achieved by reducing Bo with 
the nulling coils. A summary of reported parameters for 
various IP traps used in BEC experiments is shown in Ta- 
ble |. Typical values for the popular coil-based IP traps 
are Bo = I G, B' = 200 G/cm and B" = 150 G/cm^. 
Note that the curvature is usually much larger radially 




FIG. 2: A double-well potential with well spacing of 2zo — 
2 mm, generated by fields of Bo = - 0.75 G, B' = 200 G/cm 
and B" = 150 G/cm^ (a) The behavior for the z axis with 
the field directions indicated by arrows, (b) Field magnitude 
contours in the x-z plane with 0.1 G spacing. The barrier 
height is equivalent to a temperature {—fiB^B /kB) of 50 fiK 
for atoms in a magnetic substate with /i = /is. 



than axially so that atomic clouds are cigar shaped with 
their long axis in the z direction. For the values above, 
the radial curvature is S'.' « 40000 G/cm^ and the 
anisotropy is A = w^/wz = \JB'^ jB'!. = 16.3. In this 
calculation lOz and are the axial and radial oscillation 
frequencies in the harmonic trap. 

Ketterle et al. ^ have pointed out that a double- 
well is formed if the nulling bias (Bn) is allowed to over- 
come the pinch bias {Bp), so that Bo < in Eq. (|l|). The 
axial field now has zeros at the points z± ~ ±Zo, where 



2\Bo 



B" 



(4) 



The trapping potential then has the form of two wells, as 
shown in Fig. ^. The barrier height A_B is |i?o|, which 
occurs at (z,r) = (0,0). 

The field components, from Eq. evaluated about 
the well bottom z±, are given by 



Bx = i+B'-^B"z±) 



By - 

Bz = 



i-B' 



kB"z±) y 
B"z± z' 



(5) 



where z' 



z±. If we were to neglect the term in 



z , then these equations would be similar to those of a 
spherical quadrupole field. This approximation is rea- 
sonable near the bottom of each well, where z' /2z± <C 1. 
However, in general the curvature distorts the potential 
from a quadrupole. Note that the axial gradient scales 
with the well position Zq, so that confinement tightens 
and becomes more linear as the well separation increases. 
The gradients along the x and y axes are very similar be- 
cause we will consider well separations small enough that 
B' ^ B" Zo- The x-y asymmetry is therefore small and 
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can often be neglected, but we include it in subsequent 
calculations for completeness. 

The usefulness of the potential in Fig. || for BECs 
is severely limited by Majorana spin flips at the two 
field zeros Removing this means of trap loss 

in the double-well trap above is the motivation behind 
this work. 



II. DOUBLE-TOP TRAP 

To eliminate loss from the two quadrupole-like wells, 
we apply the TOP trap scheme developed by Cornell 
et al. [|40| for the simple quadrupole trap. We now con- 
sider adding a rotating bias field to the double-well IP 
trap. The rotating field component, with magnitude Bt, 
is chosen to be 



Bx — Bt cosLut, 
By — Bt sinwt, 
B, = 0. 



(6) 



The oscillating field has no z component, which avoids 
the introduction of further radial asymmetry and mod- 
ulation of the well spacing, stiffness, and barrier height 
during the bias rotation. This choice also leads to a rela- 
tively simple description of the trap. As in the standard 
TOP trap, the frequency of rotation lo must satisfy 



(7) 



where is the Lamor frequency in field Bt- The lower 
bound ensures a time-averaged potential, and the upper 
bound allows the atomic magnetic dipole to follow the 
oscillating field so that the atoms remain trapped. These 
limits cause no practical difficulties as trap oscillation 
and Lamor frequencies are normally of order 100 Hz and 
1 MHz respectively. 

The key mechanism of the TOP trap is that the field 
zero of a spherical quadrupole is displaced outside the 
atomic cloud by the rotating bias field of Eq. (^. In 
the original implementation pQ ] the field rotated about 
the symmetry axis of a spherical quadrupole giving a 
circular locus of i? = at radius ro — Bt/ B' , where B' is 
the radial gradient of the quadrupole. The double-TOP 
trap described in this work does not in general have equal 
gradients in the plane in which the bias rotates, so the 
locus is slightly elliptical. The axial displacements for x 
and y are given by 



Ti = 



Bt_ 

IB' I 



(8) 



where z = {x, y} and B[ ~ is the magnetic field gra- 
dient on axis i evaluated from Eq^ 
mentations of the TOP trap Q |4| 
locus because the oscillating magnetic field does not ro- 
tate about the symmetry axis of a spherical quadrupole. 
The double-TOP will have a circular locus for the field 



5^) . Some other imple- 
also have an elliptical 




(b) 50 




FIG. 3: A double-TOP potential with the same parameters as 
Fig. I {2zo = 2 mm , Bo = -0.75 G) and Bt = 1 G. (a) The 
behavior on the z-axis with numerical (solid), Mexican-hat 
(dash-dot) and double-well IP (dotted) potentials shown, (b) 
A contour plot in the x-z plane for the numerical integration 
with 0.1 G spacing starting from the well bottom Bt- The 
barrier height is equivalent to 17 fjK for atoms with jj, = fiB- 



zeros when the x-y asymmetry is negligible, correspond- 
ing to the limit Zq ^ 1 cm for our typical parameters. 

Following Cornell's method for the standard TOP trap, 
the field distribution for the double-TOP can be evalu- 
ated by averaging the field magnitude: 



ZTT 



B{t) dt, 



(9) 



where B{t) is the instantaneous field magnitude. In Fig. ^ 
we have numerically evaluated Eq. |^. The addition of the 
rotating bias field to the double- well IP trap has displaced 
the field zero as required, leaving a minimum of field 
magnitude of Bt occurring at (z, r) = {z±,0). The field at 
the top of the barrier separating the two wells is now the 
quadrature sum of static and rotating bias fields, given 
by 



Bh 



B^ + B?. 



(10) 



The barrier height AB is then given by AB ^ Bb — Bt- 
The well separation is not affected by the addition of the 
rotating bias field. 

There is no straightforward analytical solution to 
Eq. (|^) that applies for a wide range of parameters. How- 
ever, it is possible to find results for restricted conditions 
and we now present two such cases. Along the central 
axis of the trap, the field magnitude can be approximated 
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by a 1-dimensional Mexican-hat functional form 

2 

Bav(z) =Bt + AS ( 1- ( ^ 




(11) 



and the error is less than 1.5% for \Bo\/Bt < 1 in the 
region \z\/zo < V^- This description becomes more ac- 
curate as the ratio \Bo\/Bt becomes smaller. Figure ^ 
shows that the fit is best in the central region where 
Bt < Bay{z) < Bb. At larger values of z the time- 
averaged field is dominated by the static axial compo- 
nent, so the potential becomes harmonic. 

The second case is to consider the shape of each well 
near the minimum. This description applies for a cold 
atomic cloud confined within either well. We proceed 
by performing an analytical integration of Eq. (|9|), using 
the double-well fields expanded about z — ±Zo and the 
rotating bias. The resulting field magnitude is given by 
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\ z 

8Bt 



(12) 



Examining this result we find that the radial dependence 
is harmonic and is dominated by the term B'"^ / {ABt), 
which is identical to the standard TOP trap with field 
rotating in the x-y plane. The term in [y^ — x^) is the 
result of the slight x-y asymmetry of the double-well and 
is negligible for Zq <C 1 cm. Compared to the single-well 
IP trap of Eq. (H) the axial curvature has been modi- 
fied by the factor 2\Bo\/Bt, so that axial confinement is 
tightened for \Bo\ ~ Bt- The terms in z'^ and z''* de- 
scribe the barrier, and for our typical parameters these 
terms are important when Zq is smaller than approxi- 
mately 500 /^m. If the wells are more widely spaced, the 
third order term adds only a slight tilt to the harmonic 
confinement and the 4th order term can be neglected, as 
shown in Fig. 4. 

In summary, the complete potential can be calculated 
numerically, while near the well minima can be described 
analytically. The Mexican-hat functional form is a useful 
approximation for the axial field magnitude in the central 
region. 



III. LOADING FROM A lOFFE-PRITCHARD 
TRAP 

The double-TOP trap may be loaded with a Bose con- 
densate by transfer from a IP trap. To do this without 
spin-flip loss it is necessary to apply the rotating bias be- 
fore reverse biasing to form a double-well. We model this 
loading scheme by adding the rotating bias of Eq. (^) to 
the IP field components of Eq. (||) , and on integrating for 
the magnitude find that 



Sav — Bh 



1 


Bo 


(B'^ 


B"\ 


BfB'^' 


2 


B'b 


KB^ 




2B! \ 



(a) 60 

I 40 

I 20 
m 





\ 












^^^^^ 





800 



900 



1000 



1100 



1200 



(b) 30 




300 400 
z ()im) 



700 



FIG. 4: Right-hand well of double-TOP potentials with com- 
mon field parameters B' = 200 G/cm, B" = 150 G/cm^ and 
Bt = 1 G. Both the numerical average (solid) and evaluation 
of Eq. |l^ (dashed) are shown, (a) The axial potential for 
widely separated wells (zo = 1000 /^m) and relatively high bar- 
rier (A_B — 250 mG). Eq. |l^ is evaluated to second order, (b) 
The axial potential for closely separated wells (zo = 400pm) 
and a low barrier (AS = 7.2 mG). The analytical result is 
taken to 4th order in z. 
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(13) 



Compared to the initial IP trap, both the axial and ra- 
dial curvatures are reduced by at least a factor of Bo/B^ 
due to the time-average. This effect can be minimized if 
Bb ~ Bo which occurs when Bt is small in comparison 
to Bo- Since Bo > 0, there are no zeros of the field in 
this configuration and therefore there is no restriction on 
the minimum size of Bt- This choice of parameters also 
minimizes the effect of the third term in the radial cur- 
vature. A fourth order term in z is required because the 
axial curvature goes to zero with Bo- 

The loading process is an adiabatic evolution from IP 
to the time-averaged IP trap and finally to double-TOP. 
We can use Eq^ ( p^ ) to describe most of this process, as 
shown in Fig. Q The first step is to increase Bt from to 
1 G while holding Bo = 1 G. The trap minimum becomes 
Bb = ^/2Bo, and the axial and radial curvatures are re- 
duced by factors of 0.71 and 0.53 respectively. A typical 
condensate is sufficiently small (sa 150 /im long) for the 
axial confinement to be harmonic. The second stage is to 
ramp Bo down to —0.75 G which creates a double-TOP 
trap characterized by AB — 250 mG, 2zo — 2000 fim. 
During this ramp the radial confinement is mostly un- 
changed because the rotating bias dominates Bt,. When 
Bo passes through zero, Bt must be sufficiently large to 
place the field zero outside the cloud of atoms and pre- 
vent spin-fiip loss. For the typical parameters used, the 
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FIG. 5: Numerical (solid) and analytical (dash-dot) modeling 
of the transformation from IP to double-TOP trap by reduc- 
ing Bo- The behavior is shown for the z-axis (left column), 
and a;-direction at the trap minima (right column) indicated 
by the marker on the corresponding z-axis plot. The plots 
show: (a) The initial IP trap with = 1 G, B' = 200 G/cm, 
B" = 150 G/cm^ St = G; (b) The time-averaged IP formed 
with St = 1 G; (c)-(f) Bo equal to 0.25 G, G,-0.19 G, and 
—0.75 G respectively (all other parameters held constant). 
The well separations in (e) and (f) are 2zo — lOOO^m and 
2000/im. The radii of the locus traced by the two quadrupole 
field zeros in (d)-(f) are constant and equal to = 51.9 fim 
and Ty = 48.2 fim. 



radial curvature of 20000 G/cm^ at Bo = means a typ- 
ical alkali condensate of 10^ atoms will be significantly 
inside ro = 50 /im. The ramp of Bq produces a dramatic 
change in potential on the z axis. While Bo > 0, the 
curvature is positive but reducing, leading to a flattening 
as the quartic term becomes more important. Some con- 
trol of the rate of change of the curvature is possible by 
the choice of Bt relative to Bo- In this case, where they 
are approximately equal at start of the ramp, the cur- 
vature does not greatly deviate from a linear reduction, 
and this is also true for the case Bt > Bo- In the oppo- 
site extreme the curvature can collapse quite suddenly. 
When Bo passes through zero, the sign of the curvature 
reverses and the barrier starts to rise. In Fie^|(f) we see 
that although the analytical solution (Eq. (^3|)) describes 
the qualitative shape of the potential, there is no longer 
good quantitative agreement with the numerical result. 
The Mexican hat description is more appropriate for this 
situation. 

A second method for loading Bose condensates into 
a double-TOP configuration is to evaporatively cool a 
thermal cloud. In this case the evaporation process is 
much the same as for a standard TOP trap, with the 
RF field oriented perpendicular to the rotating bias field, 
except that RF resonance occurs at two points at the 
perimeter of the cloud rather than one. Alternatively, 
an atomic cloud could be cooled first in a IP trap and 
then in a double-TOP. However in this case, the standard 
orientation for the RF antenna in a IP trap is not well 
matched for optimum evaporation in the double-TOP. 



IV. CONTROLLING THE BARRIER HEIGHT 
AND WELL SPACING 

Careful control of the barrier height will be important 
in any investigation of condensate splitting and Joseph- 
son tunneling. Such studies are likely to require a bar- 
rier height comparable with the condensate chemical po- 
tential and therefore much smaller than that which we 
considered in the sections above. Recall t hat the bar- 
rier height is given by AB ~ \/ B^ + Bf — Bt, so for 
the trap parameters associated with Fig. ||(f) we have 
Ai? = 0.25 G. For atoms trapped in states with a mag- 
netic dipole moment of /^b, this barrier height is equiva- 
lent to a temperature of fisAB /ks = 16.8 fiK- In com- 
parison, the chemical potential, /i [Q, of 10^ atoms of 
^^Rb in the \F — 2,mF — "2) state and confined in one of 
the (almost harmonic) wells is /i/Zcs = 180 nK, so that 
the ratio is fiB^B/fi — 93. 

The options for lowering the barrier height are to lower 
Bo or increase Bt . Increasing Bt keeps the well positions 
fixed but relaxes the radial confinement, while decreasing 
\Bo\ reduces the well spacing. In Fig. ||(a) we show the 
change in barrier height as a function of Bo while keeping 
Bt fixed. Initially the barrier is some 100 times higher 
than the chemical potential but decreases to a compara- 
ble height at Bo = —55 mG. The chemical potential is 
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Bj(G) 

FIG. 6: (a) Control of the barrier height by ahering the static 
bias field Bo with = 1 G. In the range shown, the ratio 
of barrier height and chemical potential drops from 100 to 1. 
The circle indicates the parameters used for the inset. Inset: 
the magnetic potential for Bo — —60 mG and the chemical 
potential (horizontal lines) for 10® Rb atoms in each well. The 
calculation of fi is approximated by neglecting any efi'ect of the 
condensate in the other well, (b) Control of the barrier height 
by altering the rotating bias field Bt with Bo — —0.75 G. 
Note that even at 25 G the barrier height is still considerably 
greater than the chemical potential. 



estimated by a numerical method where we integrate over 
the magnetic potential in the Thomas-Fermi approxima- 
tion iQ. We keep the number of atoms constant at 10^ 
for each value of Bo and solve for the chemical potential. 
Note that achieving very low barriers with this method 
requires very small values of Bo, which is discussed below. 

We could also lower the barrier height by increasing 
Bt, as shown in Fig. ^(b). For our choice of parameters, 
using this method to make the barrier height equivalent 
to the chemical potential would involve increasing Bt to 
approximately 1 x 10'^ G, and is therefore impractical. 
However, with a lower value of Bo, increasing Bt offers 
greater sensitivity. For example, with Bo = —110 mG the 
barrier height and chemical potential become comparable 
when Bt is approximately 20 G. 

Achieving small well spacing also requires a very 
small bias field. For example, to obtain a well spac- 
ing of 2zo = 300 /zm requires a static bias field of only 



-17 m G for our choice of axial curvature. Note that 
Zo — \/2\Bo\l B" , so that best control is achieved for 
largest possible curvatures and widely separated wells. 

Achieving tight control of close well spacing and low 
barrier height therefore requires a high degree of coil cur- 
rent stability and shielding from external fields. High 
current power supplies of the sort often used for IP traps 
have relatively poor current stability, so that achieving 
close well separation with these is likely to be difficult. 
However, a low current example where a high degree of 
stability has been achieved is the QUIC trap ||2^ . The au- 
thors report in Ref. reducing "residual fluctuations 
in the magnetic field to a level below 0.1 mG". Since 
the bias field required to achieve control over the well 
separation is small it should be possible to add a small 
external bias field without compromising the stability of 
a QUIC-style trap. It is worth noting that in the limit 
of a very small well spacing it becomes increasingly dif- 
ficult to form a significant barrier, because this requires 
Bt ^ Bo- This condition leads to a very small Tq and 
subsequent trap loss from Majorana spin-flips. Ensher 
has also found evidence that an atomic cloud in a 
TOP trap with small Bt is vulnerable to transitions to 
untrapped magnetic spin states, induced by residual AC 
magnetic fields associated with current noise. 



V. CONCLUSIONS 

We have developed a theoretical description for a 
double-well magnetic trap suitable for confining Bose- 
Einstein condensates and cold thermal clouds. The sep- 
aration of the wells is controlled by the static bias field 
Bo, while the barrier height also depends on the rotating 
bias field Bt . We have developed analytical forms of the 
potential that describe one or both wells, which are in 
good agreement with numerical simulations of the time- 
averaged field. This trap has the convenient feature that 
it is based solely on magnetic fields and therefore avoids 
the issue of the relative stability of magnetic and opti- 
cal fields addressed in . It also offers a barrier height 
and well spacing that are tunable, so that a wide range 
of condensate phenomena may be accessible. Wide well 
separations are readily achieved and may be useful for 
interferometric applications, where a condensate in one 
well is perturbed and a second (local oscillator) is not. 
The double-TOP scheme also has the advantage that it 
is based on existing trapping technologies. 
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